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CL|. Abstract 

. Theoretical prediction of oscillations of cumulant moments of parton 

\ multiplicity distributions inside a jet supported by experimental data in 

• some multiple production processes asks for analysis of the phenomenon 

. for the whole set of available reactions. We have found out that the oscil- 

J> I lations persist in any kind of processes and increase for particles with more 
complicated structure i.e. in the order of ee, eh, hh, hA, AA. The effect is 

^ ■ not strongly dependent on the available phase space. 

I Theoretical values of moments for quark and gluon jets up to 5th rank 

I are shown. Zeros of the truncated generating function and singularities of 

Q\ • the total generating function are discussed. 

1 Introduction 

(D : 

' Multiplicity distributions and their moments contain in the integrated form all the 

^ ■ information about a multiple production process. At the same time, their analysis 

^ . is rather simple since it does not ask for multidimensional plots and its results are 

easily interpreted. It is especially important that QCD provides quite detailed 
predictions on the behavior of parton multiplicity distributions. Even though 
the direct comparison with experimental data is impossible without Monte Carlo 
simulations with definite assumptions about hadronization, qualitative features 
of the predicted (parton) and measured (hadron) distributions are so similar that 
one is tempted to say that experiment supports the predictions. Further support 
is given to the local parton-hadron duality hypothesis as well. 

The most impressive prediction of QCD is the minimum of cumulant mo- 
ments (or their ratio to factorial moments) at the rank q ^ 5 with a negative value 
at the minimum and subsequent oscillations (for review see 0). It concerns 
parton distributions. With the local parton-hadron duality at work, one could 
hope for the similar behavior of hadron distributions in gluon and quark jets i.e. 
primarily in e~^e~ annihilation. It is astonishing that such features have been 
observed not only in e~^e~ but in pp |^ , pA and AA as well. Here we complete 
the list of available reactions by analyzing HI data on ep deep inelastic processes 
P] and show regularities which appear for the whole set of the processes. 
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The similarity of patterns of cumulant moments oscillations in all the processes 
reveals some general origin of dynamics, probably indicating their underlying 
cascade nature. This similarity is especially surprising because the shapes of 
multiplicity distributions for these processes differ drastically. The minimum of 
cumulant moments at g ~ 5 becomes deeper for colliding particles with more 
complicated structure i.e. in the order of them like e^e~ ,eh, hh, hA, AA. This 
regularity has no physical interpretation yet. It does not contradict, however, to 
results obtained by Monte Carlo simulation schemes. 

In view of the existing and forthcoming ^ attempts to study multiplicity 
distributions inside gluon and quark jets separately we provide our predictions 
for their cumulant moments up to the rank 5, comparing them to DELPHI data 



We have also calculated the location of zeros of the truncated generating 
function for ep collisions and compared it with corresponding pattern in other 
reactions. The possible scaling of zeros location and the nature of the singularity 
of the total generating function are discussed. This is appealing in view of the 
fact that asymptotically the singularity is located in the vicinity of the point 
where moments are calculated as derivatives of the generating function. The 
similarity to the partition function properties in statistical machanics where they 
are related to phase transition problems gives further impact to that study. 

2 Moments of multiplicity distributions 

We have analyzed new data of HI Collaboration [^] on multiplicity distributions in 
deep inelastic scattering for different energy (W) intervals in terms of moments. 
This method proved to be useful for ee,pp 0] and hA,AA processes. The 
cumulant moments are especially sensitive to tiny details of distributions and 
reveal the intriguing oscillations. 

Let us define the normalized moments of the multiplicity distribution Pn{y) 
where P„ denotes the probability of n-particle production process at c.m.s. energy 
s^/^ and y = ln(s^/^/m) with m ~1 GeV. After measuring Pn{y) one can calculate 
the factorial moments Fq of the rank q as 



and subsequently the cumulant moments Kg by iterative procedure from the set 
of relations 
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are binomial coefficients. For all high energy experimental distributions facto- 
rial moments are larger than 1 but cumulant moments oscillate. Usually, both 
factorial and cumulant (in modulo) moments increase fast with their rank q. 
Therefore, it is convenient to consider their ratio 

= (4) 

It is especially so because QCD provides results just in terms of Hg predicting 
specific oscillations as a function of q with first minimum located at 

^ +i^5 (5) 



hi'jo 2 

where hi = 11/24 is given by the integral over the kernel of QCD evolution 
equations, and 7o = 6as/n,as is the QCD "fine structure constant" (70 ~ 0.48 
at Z^-peak). This estimate is obtained for parton distributions in gluon jets. For 
quark jets it could shift to g ~ 7. 

In Fig.l we have plotted the dependence of ratios Hg on their rank q for 
various initial colliding partners (ee, ep,pp,pA, AA) at high energies. All of them 
reveal the oscillations with first minimum located at g ~ 5 or nearby as predicted 
by QCD (for parton distributions in gluon jets!). The mentioned regularity in 
the depth of first minimum is clearly seen. It becomes deeper for those reactions 
where the internal structure of the colliding objects is more complicated. This 
regularity has been noticed in |^ and new data on ep collisions P fits it very 
well. 

The oscillations could be attributed to some common cascade nature of the 
processes even though the present theoretical approach to ee differs drastically 
from AA description (s-channel cascades versus t-channel ones). Their periodicity 
and amplitude are even more intriguing. Anyway, no conventional phenomeno- 
logical (probability theory) distribution fits them. 

The attempts to ascribe all the oscillations to the cut-off of multiplicities at 
finite energy fail because it gives rise to comparatively small values of amplitudes 
if maximal multiplicity is much larger than the average one[^] and, even more 
important, predicts smaller amplitudes in A A as compared to other processes due 
to higher multiplicities in AA. It contradicts Fig.l. Nevertheless various cut-offs 
are rather important provided they limit substantially the available multiplicity. 
The HI data at the same energy W but for different rapidity windows confirm 
this statement if plotted in terms of Hg in Fig. 2. The smaller window one has 
chosen, the larger is the amplitude of the oscillations with the first minimum 
shifted to higher ranks. 

As seen from eq. (5), the minimum location in QCD depends on coupling 
constant and on the kernel of evolution equations where the vector nature of 
gluons is of primary importance. 
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Its dependence on the QCD running coupling is such that qmin should be 
slightly larger at higher energies because 70 (or a^) decreases. This shift is rather 
hard to notice since the decrease of 70 is not strong enough and, besides, the values 
of moments are only available at integer ranks. Some qualitative indications can, 
however, be obtained if one deals with data of the same collaboration. Otherwise, 
different selection criteria are difficult to disentangle. This unique possibility is 
provided by recent HI data f^. Our analysis shows that QCD trends are rather 
weak and multiplicity cut-off effect conceals these trends. In Fig. 3 we plot the 
moments for different "excitation energies" (W) of the proton calculated at the 
same maximal multiplicity equal to 18. We see that the minimum tends to move 
to the left at higher W. Therefore, cut-off effect is still strong enough if maximal 
multiplicities are not much larger than the average multiplicity. 

Another chance to check QCD results appears when gluon and quark jets are 
accurately isolated and considered separately as is done in refs. [7| , for example. 
Then the widths and higher moments (apart average multiplicities considered in 
[Q) of multiplicity distributions within jets can be treated. In Table I, we show 
the values of corresponding moments as given by fixed coupling QCD equations. 
However the value of the coupling constant has been chosen higher than at Z^- 
peak to take into account its running during the cascade evolution. According 
to the analytic results the cumulant moments depend rather strongly on 
the coupling constant. It has been confirmed by varying its value in numerical 
estimates. The moments decrease for both gluon and quark jets if coupling 
constant increases. However we have found that the ratios H^/H^ are rather 
stable for different 70 being approximately equal to 2 for q = 2 and ranging from 
4 to 5 for q = 3. One must remember however that there are two predicted 
numbers to be compared with experiment for each rank, not a single one. 

As seen from Table I, the width of the multiplicity distribution for a quark jet 
is larger than that for a gluon jet. At present, there is DELPHI data where the 
values of the Negative Binomial Distribution parameter k for multiplicity distri- 
butions in quark and gluon jets are shown in Tables 5 and 7 of [^. Misleadingly, 
they are called dispersions in the titles to the Tables even though the higher val- 
ues of k correspond to lower dispersions. They have been derived from the fit to 
the Negative Binomial Distribution using the DURHAM and JADE algorithms. 
The values of k obtained according to JADE (Table 7 in 0) and DURHAM (Ta- 
ble 5 in 1^) algorithms differ drastically, probably due to the different angular 
coverage for those algorithms, and contradict to each other (especially, at low 
jet energies). Unfortunately, the statistical errors are shown only in the Tables. 
They are large but one can suspect that the systematical errors (not shown) can 
be even larger due to b and c quark jets not separated. To compare theory to 
the data such a separation should be done. Otherwise from the data considered 
boldly, one would conclude that there is a contradiction between theory and ex- 
periment. Better understanding of these systematic effects is needed in order to 
make a fair comparison between data and theory in this subject. Higher moments 
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ask for better jet separation and statistics also. 



3 The generating function 

Factorial and cumulant moments are calculated from experimental multiplicity 
distributions according to eqs. (l)-(3). For theorists, it is more convenient to 
deal with the generating function 

oo oo rp 

G{z, y) = E(i + ^TPniy) = E -TzHny, (6) 

where (n) = I]J^Lo ^-^n the average multiplicity. If G{z, y) is known, the facto- 
rial and cumulant moments are given by its derivatives 



^ 1 d'^G{z,y) 
^ {nY dzi 

^ 1 dnY,G{z,y) 
{n)i dz" 



(7) 



(8) 



Therefore it is important to know how close the singularity of G{z,y) in z-plane 
is located to the point z=0 where the moments are calculated. The nature of 
the singularity is also important to compare it to the generating functions of 
phenomenological distributions and to look for some statistical analogies with 
their speculations about phase transitions (see The analytic structure of 

the generating function in gluodynamics is determined by the integro-differential 
equation 

^^lpyl= rdx[--(l-x)(2-x(l-a;))]7olG(^,y+lnx)G(;^,y+ln(l-a;))-G(^,y)]. 
ay Jo X 

(9) 

In QCD it is replaced by a set of two equations. Until now, nobody was able 
to get the exact solution of this equation in terms of G{z,y). It has been solved 
in terms of the moments exactly for fixed coupling constant and in higher order 
approximations for running coupling constant (see 0]). It is well known [|T0[ that 
in the double-logarithmic (DLA) approximation the equation reduces to 



d^ \nG{z,y) 



7o'[G(z,y)-l]. (10) 



dy"^ 

It is quite clear that the leading singularity governed by eq. (10) can be written 
as 

G{z,y)^{z-zo{y))-' (11) 

because differentiating In G twice produces the second order pole in the left hand 
side that asks for the same singularity in the right hand side. Since KNO-scaling is 
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valid in this approximation, the factorial moments do not depend asymptotically 
on and the generating function (6) depends on the product of z{n) i.e. zq{ii) = 
K/{n). The constant k is equal to 2.552 |[T0[ . 



The eq. (11) reminds the negative binomial distribution with its parameter 
kNBD = 2 that corresponds to the extremely wide multiplicity distribution. In 
the last case, the constant n would, however, be equal to 2 instead of 2.552. More 
complete formula for G{z, y) near the singularity looks in DLA like 



G{z, y) = + ^ - I In ^ + 0(1). (12) 

[z-zoy Z-Zq 3 Zo 



This formula has rather heuristic than practical interest because experimental 
distributions are far from being fitted by the negative binomial distribution with 
^NBD = 2. It shows, however, that 2; = is still quite far from the singularity 
(which moves to 2; = asymptotically!) in the sense that the pole contributions 
cancel at this point in G{z, y) itself. There is no cancellation in its derivatives 
which increase like factorials due to the leading singularity. 

The attempts to go beyond DLA are connected with Taylor series expansion 
of G (see [0) so that the equation (10) gets in the right hand side additional terms 
proportional to derivatives of G{z,y). Each derivative has stronger singularity 
than the function itself and, therefore, such an approach can not be used near 
the singular points. It demonstrates also why DLA does not work in practice. 
One has to come back to the original equation (9) which has not been solved yet. 

From the experimental side, there is another way to approach the problem. 
One can study the behavior of zeros of the truncated generating function Gtr{z, y) 
which is defined by 

no 

GUz,y) = T.i^ + zrP,iy) (13) 

n=0 

and differs from the total generating function by the sum being truncated at 
some value oi n = no that happens inevitably in experiment due to finite energy. 
In principle, one can truncate at different values no for a given energy y up to 
the maximum available multiplicity Umax- The truncated generating function is a 
polynomial of order no with uq complex conjugate zeros in the complex z-plane. 
It has no singularity at finite z. The positions of zeros have been studied in 
several papers (e.g. |^, |Tl|, |12[)- It has been shown that at large uq their pattern 



reminds the circle. We have analyzed DESY data and have found that the similar 
regularity appears there as well with some additional substructure in the left hand 
side of the circle due to low values of cut-off uq. It can be seen in Fig. 4. One can 
compare our findings with previously pubhshed results [^, 0, |T^ . 



With no increasing, zeros tend to the real axis and their limit at no — > C)0 
should coincide with the singularity of the generating function. The relation of 
the density of zeros and of their approach to the real axis to the singularity nature 



has been discussed in |jT3[. This kind of analysis of experimental results has not 
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been done yet. It requires special study and we hope to do it later and publish 
elsewhere. 



4 Conclusions 

In previous papers @, ^, |ll], |12] the analysis of the moments of multiplicity 



distributions and of zeros of the truncated generating functions has been done 
in ee,pp,pA, AA collisions. We complete this list of reactions doing the same for 
deep inelastic ep coUisions using recently published HI data. We find out that this 
reaction fits quite well the general regularities observed previously for oscillations 
of the moments and zeros of the truncated generating function. In particular, 
the depth of the first minimum of oscillations is intermediate between ee and pp 
minima supporting earlier conclusion about its dependence on the structure of 
the colliding objects. 

Some predictions for quark and gluon jets are discussed in connection with 
recent experimental data and forthcoming attempts to measure jet moments. 

Also, we discuss and try to attract attention to the problem of singularities 
of the generating function which is still waiting for its solution. Our findings are 
encouraging in that respect. 
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Figure Captions 

Fig.l The ratio Hq for e^e~ , ep, pp, hA and AA collisions. (See text for more 
details) . 

Fig.2 The ratio Hq for ep collisions at 150 < W < 185 GeV for different pseu- 
dorapidity intervals. 

Fig.3 The ratio Hq for ep collisions for different W and fixed maximum multi- 
plicity equal to 18. 

Fig.4 The zeros of the truncated generating function Gtr{z) for ep collisions at 
150 < < 185 GeV. 
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Table 1: r = 1.8; np = 4; 70 = 0.648; 7 = 0.476. 
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